arXiv:astro-ph/9701245vl 30 Jan 1997 


The Density and Temperature of Molecular Clouds in M33 


Christine D. Wilson^ 

Constance E. Walker^ 

Michele D. Thornley^ 

Submitted to the Astrophysical Journal: July 23, 1996 
Revised: January 21, 1996 


Received 


accepted 


^Department of Physics and Astronomy, McMaster University, Hamilton, Ontario L8S 
4M1 Canada 

^Department of Astronomy, University of Arizona, Tucson AZ 85721 U.S.A. 
^Department of Astronomy, University of Maryland, College Park, MD 20742 U.S.A. 





2 


ABSTRACT 

We have observed the J=2-l, J=2-l, and J=3-2 lines in a 

sample of seven giant molecular clouds in the Local Group spiral galaxy M33 
using the James Clerk Maxwell Telescope. The ^^CO/^^CO J=2-l line ratio 
is constant across the entire sample, while the observed J=3-2/J=2-l 

line ratio has a weak dependence on the star formation environment of the 
cloud, with large changes in the line ratio seen only for clouds in the immediate 
vicinity of an extremely luminous HIT region. A large velocity gradient analysis 
indicates that clouds without HII regions have temperatures of 10-20 K, clouds 
with HII regions have temperatures of 15-100 K, and the cloud in the giant HII 
region has a temperature of at least 100 K. Interestingly, the giant HII region 
appears capable of raising the kinetic temperature of the molecular gas only for 
clouds that are quite nearby (< 100 pc). The continuous change of physical 
conditions across the observed range of star formation environments suggests 
that the unusual physical conditions in the cloud in the giant HII region are 
due to post-star formation changes in the molecular gas, rather than intrinsic 
properties of the gas related to the formation of the giant HII region. The 
results from this study of M33 suggest that similar observations of ensembles 
of giant molecular clouds in more distant normal spiral galaxies are likely to 
give meaningful measurements of the average physical conditions inside the 
molecular clouds. These results also imply that clouds with a factor of three 
difference in metallicity have similar density and temperature, which in turn 
imply that the differences in the CO-to-H 2 conversion factor seen in these clouds 
can be attributed entirely to metallicity effects. 


Subject headings: HII regions - galaxies: individual (M33) - galaxies: ISM 



galaxies: Local Group - ISM: individual (NGC 604) - ISM: molecules - stars: 
formation 
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1. Introduction 

Determining the physical conditions inside molecular clouds is important for 
understanding the link between the properties of the molecular gas and the types and 
amount of stars that are formed. Cloud properties that could affect the star formation 
process include the temperature and density of the molecular gas, as well as the mass 
fraction in high density gas. For example, higher gas temperatures might be required to 
form high-mass stars (Turner 1984), while a cloud with a higher mass fraction of dense 
gas might form stars with a higher efficiency (e.g. Lada et ah 1991). In return, star 
formation, particularly massive star formation, can affect conditions inside molecular clouds 
by compressing the gas at the boundaries of stellar wind or supernova shocks and heating 
the gas by increasing the ultraviolet radiation field. Giant HII regions are particularly 
interesting targets because they concentrate many hundreds of OB stars in a relatively small 
volume of space and represent the only nearby prototypes for the extreme star formation 
conditions in starburst and ultraluminous galaxies. 

In our own Galaxy, observations of low-luminosity HII regions suggest that single 
high-mass stars can form from relatively low-mass molecular clouds (Hunter & Massey 
1990). Observations of several nearby giant molecular clouds show significant differences 
in the spatial distribution of star formation. In the Taurus-Auriga region, star formation 
occurs throughout the cloud as single stars or in small clusters (Kenyon et al. 1990), while 
Orion B has most of its star formation occurring in four massive dense cores (Lada et 
al. 1991), and Orion A and Ophiuchus are each forming most of their stars in a single 
dense cluster (McGaughrean & Stauffer 1994, Wilking, Lada, & Young 1989). Despite the 
differences in the spatial distribution of the young stars, the star formation efficiency is 
1-4% for all four regions (Evans & Lada 1991). Interestingly, the star formation efficiency 
in two giant HII regions in M33 is also in the range of 2-5% (Wilson & Matthews 1995). 
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The clumpy structure of molecular clouds (e.g. Stutzki & Giisten 1990) implies 
that the volume-averaged density of a cloud is not a good measure of the density of the 
molecular gas producing most of the observed emission. For extragalactic molecular clouds, 
their low and usually unknown hlling factor within the beam also implies that the kinetic 
temperature cannot be estimated from the peak brightness temperature of the observed 
line. Instead, the density and temperature must be determined using radiative transfer 
models and observations of many line ratios of CO and its isotopomers. Because this 
technique requires observations of some of the rarer isotopomers such as ^^CO, it was hrst 
applied to starburst galaxies with strong CO lines (e.g. Tilanus et al. 1991, Giisten et 
ah 1993). However, one difficulty with even the nearest starburst galaxies is that the 
observations measure the emission from many molecular clouds within a single beam. For 
example, at a distance of 3 Mpc, a 20" beam subtends a diameter of 290 pc, which is much 
larger than the 10-100 pc diameter of giant molecular clouds in the Milky Way (Sanders, 
Scoville, & Solomon 1985). If the physical conditions in molecular clouds vary signihcantly 
from one cloud to the next, the results of htting the average line emission may bear little 
resemblance to even the average properties of the molecular cloud ensemble. 

Local Group galaxies offer an advantage for such studies because they are close 
enough (50 kpc - 1 Mpc) that it is possible to isolate individual molecular clouds within a 
15-20" beam. In addition, individual clouds can cover a substantial fraction of the beam 
and thus the GO lines are reasonably strong. Observations of a large sample of individual 
clouds can identify variations in the density and temperature from one cloud to the next, 
as well as correlations in cloud properties with the presence and intensity of massive star 
formation. The uniformity (or lack thereof) of the molecular cloud population also provides 
a means to test the reliability of radiative transfer techniques for determining density and 
temperature in more distant galaxies, where only the average emission from the molecular 
cloud population can be observed. In addition, density and temperature measurements 



6 


for individual clouds can be used to assess whether changes in the density and brightness 
temperature of the gas are likely to produce systematic errors in the observed correlation of 
the CO-to-H 2 conversion factor with metallicity (Wilson 1995, Dickman, Snell, & Schloerb 
1986, Sakamoto 1996). 

In this paper we present observations of seven giant molecular clouds in the spiral 
galaxy M33. The clouds were chosen to cover a wide variety of star formation conditions, 
from clouds with no optical HII regions to a cloud located in the brightest giant HII region 
in the galaxy. The clouds also cover a range of three in oxygen abundance (Vilchez et ah 
1988). The observations and data reduction are discussed in §2. The observed line ratios 
are compared with previous Galactic and extragalactic observations and with the star 
formation environment of the clouds in §3. The density and temperature obtained from 
an analysis of the line ratios with a large velocity gradient code are presented in §4. The 
implications of the results for observations of more distant galaxies and the calibration of 
the CO-to-H 2 conversion factor as a function of metallicity are discussed in §5. The paper 
is summarized in 56. 


2. Observations and Data Reduction 

Seven giant molecular clouds in M33 were observed with the James Clerk Maxwell 
Telescope (JCMT) in the J=2-l, J=2-l, and J=3-2 rotational transitions 

over four separate observing runs (1993 July 30 - August 1; 1993 August 25-26; 1994 July 
16-18; and 1995 July 31 - August 1). The half-power beamwidth of the JCMT is 22" at 
230 GHz and 15" at 345 GHz. The clouds were selected from the interferometric samples 
of Wilson & Scoville (1990, 1992) and have also been observed in the ^^CO and ^^CO 
J=l-0 lines with 55" resolution (Wilson & Walker 1994). One position (MG 1) includes 
emission from two other clouds near the half-power point of the 22" beam (MG 2, MG 
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4, Wilson & Scoville 1990), while only a single cloud produces the CO emission at the 
other positions. Each cloud was observed at a single position in the J=2-l and 

J=2-l lines and in a hve-point cross spaced by 8" in the J=3-2 line. All observations 

were obtained by position-switching to a location in the inner galaxy (a(1950)=01:31:16.9, 
5(1950)=30:24:54) that had been determined previously to be free of J=l-0 and J=2-l 
emission. Observations of the ^^CO J=2-l line were obtained at two or more different 
velocities per source to avoid the problem of spikes in the spectrometer masquerading 
as a weak line. Except for the observations in 1993 July, which were obtained with the 
Canadian Acousto-Optical Spectrometer, all observations were obtained using the Dutch 
Autocorrelation Spectrometer. A log of the observations is given in Table |l|. 

The calibration was monitored by frequently observing both planets and spectral line 
calibrators. The main beam efficiencies determined by observing Saturn, Mars, and Jupiter 
during the runs in 1993 and 1995 agreed well with the published main beam efficiencies, and 
thus for these runs we adopt the main beam efficiencies from the JCMT User’s Guide of 0.69 
at 220 and 230 GHz and 0.58 at 345 GHz. The spectral line calibrators were the evolved 
stars GRL 2688, GRL 618, and the compact HII region NGG 7538 IRSl. The calibrator 
observations agree very well with reference spectra available at the JGMT, with an average 
ratio of the observed to reference peak temperature of 1.00 with an rms scatter of 13%. The 
rms dispersion of the individual measurements is 10% for the ^^GO observations and 20% 
for the ^^GO observations. We attribute the larger scatter in the ^^GO observations to the 
weakness of the lines and adopt 10% as our absolute calibration uncertainty in any single 
line on the basis of the ^^GO observations. 

Throughout the 1994 run, the peak temperatures observed for GRL 618 and GRL 
2688 were consistently 70% of the values measured in the other runs. The low efficiency 
of the telescope was conhrmed by remeasuring the ^^GO J=2-l spectrum of MG 20, which 


was also 70% of the previously measured value. The main beam efficiency determined 
from observations of Mars and Uranus was also low by the same amount. Thus for the 
data obtained during this run, we adopt a main beam efficiency at 220 and 230 GHz of 
0.70 X 0.69 = 0.48. 

As we are comparing line strengths measured with the same beam diameter (either 
the original beam or convolved to a larger beam), the best temperature scale to use would 
be the scale. Although the clouds in M33 are much smaller than the full diffraction 
pattern of the telescope (5-16"), the scale corrects for all telescope efficiencies except 
the source-beam coupling efficiency, which should be the same at all frequencies since we 
have observed the emission with the same beam at each frequency. However, to convert 
the observed temperatures to requires knowing rjfgs, the forward scattering and 
spillover efficiency. This efficiency is difficult to measure, and so we did not attempt to 
measure r]fss during our observing runs. As a result, we do not know how to correct the 
published values of r]fss for the calibration problems encountered in the 1994 run. However, 
we do have good measurements of the main beam efficiency for each run and under normal 
calibration conditions, adopting the main beam temperature scale would only change the 
J=3-2/J=2-l line ratios by 5% compared to the scale. Thus we use the main beam 
temperature scale throughout this paper, and in particular in calculating the line ratios for 
comparison with radiative transfer models. 

Given the angular size of the M33 molecular clouds (5-16"), the main beam temperature, 
Tmb, is a better approximation to the brightness or radiation temperature, Tr, than is the 
antenna temperature, T^. The three temperature scales are related by 

rp _ _ TmbVmb 

J-R — - — - 

VcVfss VcVfss 

where r]c is the source-beam coupling factor, rjfss is the forward scattering and spillover 
efficiency, and r]MB is the main beam efficiency. We can estimate r]c for the M33 clouds since 
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they have been observed at high resolution (Wilson & Scoville 1990, 1992). For the 230 GHz 
observations, rjc > 0.9, while rjfss = 0.80 and timb = 0.69, so that Tr ~ O.QTmb — l-OT^s- 
At 345 GHz, rjc > 0.8, while rjfss = 0.70 and tjmb = 0.58, so that Tr ~ 0.8Tms — I.OTms- 
We use the main beam temperature scale throughout this paper, and in particular in 
calculating the line ratios for comparison with radiative transfer models. 

The data were reduced using the Bell Labs data reduction package GOMB. The data 
were binned to a resolution of 1 km s“^ and hrst to third order baselines removed. The 
lines are quite narrow (~ 20 km s“^) and the spectrometer bandwidth was always at least 
200 km s“^, so the higher order baselines should not introduce signihcant errors into the 
hnal line intensities. Each set of hve ^^GO J=3-2 spectra were then convolved to simulate 
a 22" beam. The spectra for each cloud are shown in Figure |^. Ghannel-by-channel 
ratios of the two spectra were used to measure the average line ratio for each cloud. Only 
channels with a signal-to-noise ratio of at least 2 in the line ratio were used. This restriction 
resulted in the ^^GO/^^GO J=2-l line ratio being obtained from the central 5-14 km s“^ 
(0.6 — 1.3 X AVfwhm), while the ^^GO J=3-2/J=2-l line ratio was obtained from the 
central 9-19 km s“^ (1.1 — 1.7 x AVfwhm)- The measured line ratios and ^^GO J=2-l 
integrated intensity for each cloud are given in Table ^ 


3. CO Line Ratios and the Star Formation Environment 

The ^^GO/^^GO J=2-l line ratios given in Table ^ have an average value of 7.3 and 
an rms dispersion of 1.3 (18%). Thus we conclude that there is no intrinsic variation in 
the ^^GO/^^GO J=2-l line ratio from one cloud to another in M33. This average line 
ratio compares quite well with values measured in other galaxies such as M51 (6-10 in 
center and spiral arms, Garcia-Burillo, Guelin, & Gernicharo 1993) and NGG 4414 (9 ± 1, 
Braine, Gombes, & van Driel 1993). The value in M33 also agrees quite well with the 
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value measured in a pencil-beam survey of the Milky Way (5.5 ± 1,Sanders et al. 1993) 
and in large-area maps of the Orion A and B molecular clouds (4-6 along the main ridge, 
Sakamoto et al. 1994). The M33 value is slightly smaller than the value measured in the 
starburst galaxy M82 (12 ± 3, Tilanus et al. 1991) and in a large sample of starburst and 
interacting galaxies (13 ± 5, Aalto et al. 1995). Previous measurements have indicated that 
the ^^CO/^^CO J=l-0 line ratio is also somewhat larger in starburst and merging galaxies 
compared with other galaxies (Becker & Freudling 1991; Aalto et al. 1991; Combes et al. 
1991; Casoli et al. 1991). 

The J=3-2/J=2-l line ratios given in Table ^ show more scatter. In particular, 

the line ratio obtained for NGC 604-2 (1.07 ± 0.03 is signihcantly higher than the value for 
any other cloud. Excluding NGC 604-2, the average ^^CO J=3-2/J=2-l line ratio is 0.69 
with an rms dispersion of 0.15 (21%). This average line ratio for the six clouds in M33 is 
in good agreement with values measured in other galaxies such as the starburst galaxies 
M82 (0.8 ± 0.2, Glisten et al. 1993), NGG 253 (0.5 ±0.1 in the disk. Wall et al. 1991), and 
IG 342 (0.47 ± 0.12, Irwin & Avery 1993). The average value for the six clouds in M33 
also agrees well with values measured in M51 (0.7 in the spiral arms, Garcia-Burillo et al. 
1993) and in a pencil-beam survey of the Milky Way (0.55 ± 0.08, Sanders et al. 1993). It 
also agrees with a previous measurement for M33 in a slightly larger beam (0.64 ± 0.28, 
Thornley & Wilson 1994). Higher line ratios similar to the value for the cloud NGG 604-2 
are seen in the nuclei of both NGG 253 (~ 1, Wall et al. 1991) and M51 (1.1, Garcia-Burillo 
et al. 1993). 

The clouds observed in M33 have oxygen abundances that range over a factor of 
three (Vilchez et al. 1988). If we momentarily ignore the unusually high ^^GO J=3-2/2-l 
line ratio in NGG 604-2, there is no evidence that the observed line ratios depend on the 
metallicity of the cloud. In particular, the two line ratios for MG 19 and NGG 604-4 are 
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very similar despite a factor of three difference in the oxygen abundance (Table NGC 
604-2 and NGC 604-4 are located within 120 pc of each other and presumably share the 
same oxygen abundance, so we attribute the difference in the ^^CO J=3-2/2-l line ratios 
for these two clouds to some factor other than metallicity, most likely the star formation 
environment. 

If we examine the star formation environments of our cloud sample, two of the clouds 
(MG 19 and MG 32) he between 100 and 150 pc away from the nearest optical HII region 
(Wilson & Scoville 1991), while three clouds (MG 1, MG 13, and MG 20) contain optical 
HII regions. The cloud NGC 604-2 is located near the center of the giant HII region NGC 
604, on the southern edge of the OB association which presumably powers the HII region 
(Wilson & Matthews 1995), while the cloud NGC 604-4 is located at the very edge of the 
southern extent of the Ha emission from the HII region. If we compare the CO line ratios 
for the three clouds without HII regions (MG 19, MG 32, and NGC 604-4) with those for 
the clouds with normal HII regions (MG 1, MG 13, and MG 20), the average ^^CO/^^CO 
J=2-l line ratios for the two sets of clouds agree very well. The ^^CO J=3-2/J=2-l line 
ratios for the clouds with optical HII regions (0.79 ± 0.05) are somewhat higher than those 
for the clouds without optical HII regions (0.58 ± 0.06), while the line ratio for NGC 604-2 
in the giant HII region is even higher (1.07 ± 0.03). The difference between the two clouds 
(NGC 604-2 and NGC 604-4) in and around the giant HII region is striking. The projected 
separation of the two clouds is only 120 pc (30"), yet only the cloud in the immediate 
vicinity of the giant HII region and its powering OB association has an unusual line ratio. 

The high ^^CO J=3-2/J=2-l line ratio in NGC 604-2 may be a clue to unusual physical 
conditions in the molecular clouds from which the HII region NGC 604 formed, i.e. a 
pre-star formation difference in the physical conditions that may have caused the formation 
of the giant HII region. Alternatively, the high line ratio may be due to heating of the 
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gas by the massive stars, i.e. a post-star formation change in the physical conditions in 
the molecular cloud. This interpretation implies that the giant HII region has a relatively 
small sphere of influence over which its intense radiation held can change the properties 
of the dense molecular gas. It appears that both an unusually intense radiation held and 
a cloud in close proximity to the source of the ionizing radiation are required to produce a 
large change in the CO line ratios. However, the J=3-2/J=2-l line ratios observed for 
the clouds with optical HII regions are higher than those for clouds without HII regions 
at the 2a level. This diherence suggests that HII regions that are more than an order of 
magnitude less luminous than NGC 604 can produce smaller changes in the CO line ratios. 
Thus the data for M33 provide some evidence for a continuous change of line ratio with 
increasing intensity in the star formation environment. 

One complication to this interpretation is that the optical HII regions could in principle 
be separated from the molecular cloud by up to the scale height of the disk (~ 200 pc) 
and simply appear projected on the cloud. However, the fact that most HII regions appear 
close to molecular clouds (Wilson & Scoville 1991) suggests that there is a true physical 
association between the clouds and the HII regions. Also, the three clouds without optical 
HII regions could actually be forming stars, since the clouds have sufficient optical depths 
to hide an HII region if it were located on the far side of the cloud. However, since the 
^^CO lines are optically thick for these clouds, we can only see emission from the near and 
presumably cool side of the cloud and so whether or not the clouds have an HII region 
on their far side is irrelevant to this analysis. The increasing line ratio with increasing 
HII region luminosity is consistent with a picture in which the ultraviolet radiation held 
of the optically visible HII region heats the molecular gas. We conclude that the ^^CO 
J=3-2/J=2-l line ratio of a cloud has a weak dependence on the star formation environment 
of the cloud, with large changes in the line ratio seen only for clouds in the immediate 
vicinity of an extremely luminous HII region. 
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4. Density and Temperature Determined from Large Velocity Gradient 

Models 

4.1. The Large Velocity Gradient Model 

There are several different assumptions that can be made in calculating source 
parameters from molecular line data. The most common assumption is that the gas is in 
local thermodynamic equilibrium (LTE), where the number density of molecules in any 
particular energy level is only a function of temperature. Probably the next most common 
assumption is that there is a large velocity gradient (LVG) in the gas cloud (e.g. Scoville 
& Solomon 1974). The presence of such a gradient has the effect of Doppler-shifting the 
emission of a molecule relative to that of molecules elsewhere in the cloud. When the mean 
free path of a line photon exceeds the local velocity dispersion divided by the velocity 
gradient, it can travel through the cloud without being absorbed elsewhere. A third 
approach is to assume that the velocity held is microturbulent. In microturbulent radiative 
transfer models, the velocity structure is due to thermal motions and small-scale random 
motions (White 1977). Unlike the LVG case, neighboring cloud components are radiatively 
coupled, with the degree to which they are coupled depending on the magnitude of the 
random motions and the amount of resonance scattering. A fourth possibility is to use a 
photon mean escape probability formalism to simplify radiative transfer calculations. In 
this approach, level populations and the internal intensity are calculated for an “average” 
location in a uniform, spherical cloud, with the photon escape probability from this position 
being determined by the optical depth (Black & Aalto 1991, private communication). 
These latter three non-LTE techniques often yield similar results when applied to galactic 
molecular clouds (White 1977) and thus we have chosen to use the LVG formalism in this 
analysis. 


To estimate the physical conditions in the molecular clouds in M33, we used the 



14 


large velocity gradient code RAD written by Lee Mundy and implemented as part 
of the MIRIAD data reduction package. Models were run for kinetic temperatures 
Tk = 10,15,20,30,50,100,200,300 K and for three different values of the [^^CO]/p^CO] 
abundance ratio (30, 50, 70) to span the range of abundances seen in the Milky Way 
(hanger & Penzias 1990). For each kinetic temperature and abundance combination, the 
H 2 density range was nH 2 = 10 — 10® cm“^ and the column density range per unit 

velocity was A^(^^CO)/AV = 10^® — 10^° cm“^ (km s“^)“L We use the full-width half 
maximum velocity to calculate the column density from N/AV, where the average 

velocity width for the six clouds excluding NGC 604-2 is 9.3 ± 0.6 km s“^. 

One drawback of both the LTE and non-LTE techniques described here is that they do 
not take into account the effects of cloud-cloud shielding. Cloud-cloud shielding occurs when 
two or more molecular clouds are aligned along our line of sight in both space and velocity. 
In this case the radiation from the background cloud(s) may be fully or partially absorbed 
by the foreground cloud. When a pair of isotopes, such as and ^^CO, are being used 

to probe the optical depth along a line of sight through a galaxy, cloud-cloud shielding 
will tend to produce overestimates of the optical depth in individual clouds. Optical depth 
estimates through the ensemble of clouds are only affected by cloud-cloud shielding when 
the rarer isotope {e.g. ^^CO) becomes optically thick and so does not sample all the clouds 
along the line of sight. Cloud-cloud shielding is unlikely to be a signihcant problem in M33, 
which has a relatively low molecular gas surface density and is not highly inclined. 


4.2. Model Results for M33 Molecular Clouds 

To provide better constraints for the models, we include the ^^CO/^^CO J=l-0 line 
ratio observed for these clouds in a 55" beam (Wilson & Walker 1994). (The possibility 
that the larger beam size of the J=l-0 data can produce systematic errors in our results is 
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discussed later in this section.) Although both the NGC 604 clouds were observed within a 
single beam, we can use their signihcantly different central velocities to isolate the line ratio 
for each cloud. The average ^^CO/^^CO J=l-0 line ratio for the seven clouds is 9.7 ± 0.6 
with an rms dispersion of 1.5 or 15%. We combine this average J=l-0 line ratio with the 
average J=2-l line ratio in our modeling. Since there appears to be a trend in the 
J=3-2/J=2-l line ratio with star formation environment, we used four values of this line 
ratio in our analysis: the average value for all six clouds (0.69 ± 0.06); the average value 
for clouds without HII regions (0.58 ± 0.06); the average value for clouds with HII regions 
(0.79 ± 0.05); and the value for NGC 604-2 (1.07 ± 0.16), where the total uncertainty in 
the line ratio for NGC 604-2 has been estimated at 15% to match the expected calibration 
uncertainty in the line ratio. 

Acceptable solutions for each set of three line ratios (^^CO/^^CO J=l-0, ^^CO/^^CO 
J=2-l, and ^^CO J=3-2/J=2-l) were found by visual inspection of plots of the line ratios as 
a function of density and column density at each kinetic temperature. Acceptable solutions 
were defined to be ones for which all three line ratios agreed within their la uncertainties. 
Examples of representative solutions are given in Figure § and the full range of solutions is 
given in Table The ranges of values given for temperature, density, and column density 
for each [^^CO]/[^^CO] abundance ratio are correlated in the sense that solutions with 
higher temperatures have higher column densities and lower densities. 

The results from the LVG analysis show that the kinetic temperature increases from 
clouds without HII regions (T^ = 10 — 20 K) to clouds with HII regions {Tk = 15 — 100 
K) to NGC 604-2 [Tk > 100 K). The column density also increases by at least an order of 
magnitude from clouds without HII regions to NGC 604-2. The molecular gas in all seven 
clouds is quite dense (10^ — 3 x 10^ cm“^), although the density is lower in NGC 604-2 than 
in the clouds without HII regions. Thus an increase in the ^^CO J=3-2/J=2-l line ratio 
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produces LVG solutions with higher kinetic temperatures, higher column densities, and 
lower densities. The volume-averaged densities (mass divided by volume) are 40-210 cm“^ 
for the three clouds without HII regions and ~ 200 — 400 cm“^ for the three clouds with 
HII regions (calculated from Wilson & Scoville 1990, 1992), signihcantly lower than the 
densities derived from the LVG analysis (5 x 10^ — 3 x 10^ cm“^ and 2 x 10^ — 10^ cm“^, 
respectively). The volume hlling factor of the dense gas is 0.1 — 4% for clouds without HII 
regions and ~ 2 — 20% for clonds with HII regions. For NGG 604-2 the volume-averaged 
density is somewhat higher, ~ 500 cm“^, which combined with the somewhat lower density 
from the LVG model of 10^ — 3 x 10^ cm“^ gives a hlling factor for the dense gas of 17-50%. 
These results compare well with the volume hlling factor of 10% for clumps in Orion 
(Sakamoto et al. 1994). 

Are the unusual physical conditions in NGG 604-2 due to the presence of the giant 
HII region, or do they represent pristine conditions in the molecular gas that could have 
played a role in initiating the formation of the giant HII regions? An increase in the kinetic 
temperature of the gas could be easily accomplished by the intense ultraviolet radiation held 
produced by the dense cluster of OB stars powering the giant HII region. The gas colnmn 
density could be increased throngh shock-initiated merging of two or more molecular clouds 
or through a partial collapse of the molecular cloud due to an increase in the external 
gas pressnre. However, the molecular column density could also be decreased through 
photo-dissociation of the molecular gas and thus it is difficult to predict the net ehect on 
the column density from the formation of the giant HII region. It is difficult to understand 
how the formation of the giant HII region could have produced the higher hlling factor and 
lower density of the dense clumps in the molecular cloud. One possible scenario would be 
that the formation of the dense star cluster depleted the gas in the densest clumps, leaving 
behind the lower density inter-clump gas. This gas could then be compressed to somewhat 
higher density as the cloud shrank dne to the increased external gas pressnre. In any 
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case, the continuous range of all physical conditions (kinetic temperature, column density, 
density, and hlling factor) observed from clouds without HII regions to clouds with normal 
HII regions to this cloud in a giant HII region suggests that the different physical conditions 
in NGC 604-2 are due to post-star formation changes in the molecular gas, rather than 
intrinsic conditions of the pre-star formation molecular cloud. 

The similarity of the CO line ratios for the six clouds in the normal disk despite a 
range of a factor of three in oxygen abundance might naively suggest that the [^^CO]/[^^CO] 
abundance ratio does not change signihcantly when the oxygen abundance is varied by a 
factor of three. However, this conclusion is in contrast to the Milky Way, where hanger & 
Penzias (1990) have measured a signihcant radial gradient in the [^^CO]/[^^CO] abundance 
ratio. In fact, for kinetic temperatures of 10-20 K there do exist possible LVG solutions 
with different values for the [^^CO]/[^^CO] abundance ratio but with the same density, 
kinetic temperature, and ^^CO column density which produce identical CO line ratios 
(Figure This result can be easily understood as being due to the low optical depth of 
the ^^CO emission and the high optical depth of the ^^CO emission. As long as the ^^CO 
column density is hxed, the observed line ratios will be the same for a large range in the 
[^^CO]/[^^CO] abundance ratio. Thus our data are consistent with M33 containing a similar 
radial gradient in the [^^CO]/[^^CO] abundance ratio as found in the Milky Way, without 
requiring any conspiracy between several physical parameters to produce the observed 
constant CO line ratios. 

Using the average ^^CO J=3-2/J=2-l line ratio for the six clouds (excluding NGC 
604-2) in the LVG analysis produces solutions for temperature, density, and column density 
that are intermediate between those obtained for the clouds with and without HII regions 
separately (Table |^). The mean value for each parameter derived using the average line 
ratio is within a factor of two of the average of the mean parameters derived for clouds 


with and without HII regions separately. Thus physical conditions derived by averaging 
together emission from clouds both with and without HII regions provide a reasonable 
approximation to the average physical conditions in this group of six molecular clouds. 


4.3. Possible Sources of Systematic Errors 

There are several possible complications that could introduce systematic errors into 
the LVG results. For example, the difference in the beam area used in the ^^CO/^^CO 
J=l-0 ratio compared with the other ratios could introduce systematic errors if the J=l-0 
line ratio varies as a function of the spatial resolution used (Wilson & Walker 1994). We 
can assess what effect this would have by looking at the LVG results using a lower value of 
5 ± 1, which is similar to the values of 3 to 6 seen in Galactic molecular clouds (Gordon & 
Burton 1976, Solomon, Scoville, & Sanders 1979). For simplicity we will use the results for 
the average of all six clouds in this analysis. For a [^^GO]/[^^GO] abundance ratio of 50 and 
this lower J=l-0 line ratio, the solutions permitted for the six clouds are Tk = 20 — 300 K, 
n = 60 — 2000 cm“^, and iV(^^GO) = 9 x 10^^ — 3 x 10^® cm“^. The solutions permitted for 
NGG 604-2 are = 50 — 300 K, n = 30 — 70 cm“^, and iV(^^GO) = 1 — 2 x 10^® cm“^, with 
an additional possible solution with = 300 K, n = 600 cm“^, and iV(^^GO) = 2 x 10^® 
cm“^. Much of the increased range of permitted solutions can be attributed to the larger 
uncertainty in the lower adopted line ratio (20%) compared with the observed line ratio 
(6%). On average, the effect of adopting this lower J=l-0 line ratio is to increase the kinetic 
temperature range for which solutions are available, to decrease the density by a factor 
varying from a few to more than an order of magnitude, and to increase the column density 
and its range. 

Yet another possible difficulty is that the sizes of the clouds may be different in the 
different transitions. In applying the LVG analysis we have implicitly assumed that the 
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beam filling factor is the same in each line. Recent observations of the Orion molecular 
cloud show that the cloud appears smaller in the J=2-l line than in the J=l-0 line, 

while the J=2-l line may have an even narrower extent (Sakamoto et al. 1994). If 

the J=2-l and J=3-2 emission regions have different angular extents, we would expect 
the J=3-2 emission to originate in the smaller region. In this case, the observed 
J=3-2/J=2-l line ratios would set only lower limits on the true line ratios, and so the 
true kinetic temperatures would be larger than those derived here. In this scenario, the 
large J=3-2/J=2-l line ratio in NGC 604-2 could be due to a larger ratio of the 
J=3-2 to J=2-l emitting areas compared with the other clouds. These hypotheses could be 
tested with an imaging submillimeter interferometer such as the Smithsonian Astrophysical 
Observatory’s Submillimeter Array, which is presently under construction on Manna Kea. 

Another possible problem is that more than one temperature component or a smooth 
range of temperatures may be present in some or all of the clouds in our sample (for 
example, see models by Gierens, Stutzki, & Winnewisser 1992). This is particularly true 
for NGG 604-2, which, depending on its three-dimensional orientation relative to the 
giant HII region, could be subject to intense heating on only one side. Some support 
for a two-component model is given by the brightness temperatures predicted by models 
with Tk > 100 K. The brightness temperature predicted for the ^^GO J=2-l line is > 60 
K, which when combined with the observed peak temperature of 1.19 K gives a filling 
factor of only 2%. This filling factor corresponds to a source diameter of only 3'' in a 
22" beam, significantly smaller than the 9" diameter measured in the ^^GO J=l-0 line 
(Wilson & Scoville 1992). A crude two-component model for NGG 604-2 could consist of 
one component with line ratios like the six other clouds and a hot component with ^^GO 
J=3-2/2-l ~2. This model predicts that the two components contribute about equally to 
the ^^GO J=3-2 emission, while roughly three-quarters of the ^^GO J=2-l emission arises 
in the cooler component. This cool component gives a brightness temperature of ~ 5 K 
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for Tk = 10 K, which corresponds to a source diameter of ~ 9". Additional circumstantial 
evidence in support of more than one temperature component comes from the measurement 
of the J=2-l/l-0 line ratio in MC 32 and MC 19 (Thornley & Wilson 1994). Despite 
the relatively large measurement uncertainty (0.67 ± 0.19 in a 54" beam), this line ratio 
only overlaps the solutions for the three line ratios discussed here for Tx = 15 K and 
[^^CO]/[^^CO] abundance ratios of 50-70. A more precise measurement of this line ratio 
might easily remove this area of overlap, which would again suggest that the lowest-J 
rotational transitions originate in a cooler component than the J=3-2 transition. 


5. Implications for Distant Galaxies and the CO-to-H 2 Conversion Factor 

The relative uniformity of the line ratios of the six molecular clouds observed in the 
normal disk of M33 suggests that similar observations of ensembles of molecular clouds 
in more distant galaxies are likely to produce meaningful measurements of the average 
physical conditions of the molecular gas. This conclusion is also supported by the reasonable 
agreement between physical conditions derived using the average line ratios for the six 
clouds and the average of the physical conditions derived separately for clouds with and 
without HII regions (§3). The relatively small sphere of influence (~100 pc) of the giant 
HII region NGC 604 suggests that in normal galaxies only the most intense star forming 
regions may produce signihcant changes in the molecular gas. The change in the line ratios 
is likely to be measurable only in relatively nearby galaxies (< 10 Mpc), where the warm 
molecular gas in the HII region is not diluted by emission from cooler gas included in the 
beam. It would be interesting to test whether similarly uniform line ratios are observed in 
individual clouds in the intense ultraviolet held of a starburst galaxy, but such observations 
must await the construction of an imaging submillimeter interferometer. 


The CO-to-H 2 conversion factor, a, is predicted to depend on the density and 
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brightness temperature of the molecular clouds as y/n/Ts, as well as on the metallicity 
(Dickman et ah 1986, Sakamoto 1996). Observations of Local Group galaxies, including 
several of these clouds in M33, have been used to calibrate the dependence of the conversion 
factor on metallicity (Wilson 1995). However, given the generally lower shielding by dust 
and higher ultraviolet radiation held in low-metallicity irregular galaxies, it is possible that 
molecular clouds in lower metallicity environments have systematically different brightness 
temperatures and densities. We can use our solutions for the physical conditions in the 
molecular clouds in M33 to determine the possible range in ydr/Tg and to see whether 
systematic errors in the metallicity calibration are likely to exist. 

Ignoring the possible effects of density and brightness temperature differences, Wilson 
(1995) obtained a oc (12 + logiOj. For the six molecular clouds (excluding NGC 
604-2), ^Jn/TB can vary by a factor of six for kinetic temperature solutions from 10 to 30 
K. This large range in y/n/Ts is produced by the tendency of the solutions to have lower 
densities and higher brightness temperatures for higher kinetic temperatures. As a result, it 
is theoretically possible for NGG 604-4, which has a metallicity three times smaller than the 
other five clouds, to have ^Jn/TB a factor of six smaller, which would produce a signihcantly 
steeper slope (~ —2) in the derived metallicity calibration. However, to produce this factor 
of six difference in ^/u/Tb requires NGG 604-4 to have a [^^GO]/[^^GO] abundance ratio of 
30, while the clouds in the inner disk would need an abundance ratio of 50. Thus to produce 
signihcant systematic errors in the metallicity calibration of the GO-to-H 2 conversion factor 
would require a radial [^^GO]/[^^GO] gradient in the opposite sense from that observed 
in the Milky Way (hanger & Penzias 1990). In addition, the different densities (4 x 10^ 
cm“^ vs. 2 X 10^ cm“^) and excitation temperatures ( 30 K vs. 10 K) for NGG 604-4 
and the inner disk clouds would have to “conspire” to produce the same observed GO line 
ratios despite large physical differences. We conclude that the molecular clouds in M33 
provide no evidence for a systematic error in the metallicity calibration of the GO-to-H 2 



conversion factor due to differences in the density or brightness temperature from low to 
high metallicity clouds. 


6. Conclusions 

We have observed the J=2-l, J=2-l, and J=3-2 lines in a sample of 

seven giant molecular clouds in the Local Group spiral galaxy M33. The clouds were chosen 
to cover a range of star formation conditions, from clouds without optical HII regions to 
a cloud in a giant HII region. We hnd that the ^^CO/^^CO J=2-l line ratio is constant 
across the entire sample, while the J=3-2/J=2-l line ratio is somewhat smaller in 

three clouds without optical HII regions than in three clouds with HII regions. The seventh 
cloud, located in the brightest giant HII region in the galaxy, has an even higher 
J=3-2/J=2-l line ratio. We conclude that the J=3-2/J=2-l line ratio of a cloud has 
a weak dependence on the star formation environment of the cloud, with large changes in 
the line ratio seen only for clouds in the immediate vicinity of an extremely luminous HII 
region. 

We used a large velocity gradient code to determine the density, temperature, and 
column density for the clouds. The analysis indicates that clouds without HII regions have 
temperatures in the range of 10 to 20 K and clouds with HII regions have temperatures in 
the range of 15 to 100 K, while the cloud in the giant HII region has a kinetic temperature 
of at least 100 K. We note that the giant HII region seems to have a relatively limited 
sphere of influence within which it can heat the molecular gas, since a molecular cloud 
located only 120 pc from the giant HII region shows a normal cooler temperature. The 
column density increases by about an order of magnitude from the cool clouds to the hot 
cloud, as does the volume hlling factor of the dense gas, while the density decreases by less 
than an order of magnitude. The continuity of physical properties (kinetic temperature. 
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density, column density, and filling factor) across the range of star formation environments 
in our sample suggests that the unusual physical conditions seen in the cloud in the giant 
HII region are due to post-star formation changes in the molecular gas, rather than intrinsic 
properties of the gas related to the formation of the giant HII region. 

Excluding the giant HII region, the relatively uniform line ratios observed in M33 
suggest that average physical conditions determined from similar measurements of ensembles 
of giant molecular clouds in more distant normal spiral galaxies are likely to give physically 
meaningful results. These uniform line ratios also imply that the average brightness 
temperature and density are similar in clouds ranging over a factor of three in metallicity, 
which in turn suggests that differences in the CO-to-H 2 conversion factor seen in these 
clouds can be attributed entirely to metallicity effects. 

The research of CDW is supported through a grant from the Natural Sciences and 
Engineering Research Council of Canada. MDT thanks NRC Canada for travel support. 
The JCMT is operated by the Royal Observatories on behalf of the Particle Physics and 
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Fig. 1.— J=2-l, J=2-l, and J=3-2 spectra for seven giant molecular 

clouds in M33. The spectra are binned to 1 km s“^ resolution and scaled to the main beam 
temperature scale. All spectra refer to a 22" beam; the higher resolution J=3-2 data 

have been convolved to this larger beam. The J=2-l spectrum has been scaled up by 

a factor of 3. 

Fig. 2.— (a) A typical large velocity gradient solution for the average CO line ratios for 
six giant molecular clouds in M33. The lines indicate the observed ila values for the line 
ratios of: ^^CO/^^CO J=l-0 (9.7 ± 0.6, dashed line); ^^CO/^^CO J=2-l (7.3 ± 0.5, dotted 
line); J=3-2/J=2-l (0.69 ± 0.06, solid line). Acceptable values of molecular hydrogen 
density nH 2 and column density per unit velocity N{CO )/for a kinetic temperature 
Tk = 20 are found where all three sets of lines intersect. This hgure shows solutions for an 
adopted [^^CO]/[^^CO] abundance ratio of 50; slightly warmer (30 K) solutions are possible 
if we adopt an abundance ratio of 30. (b) A similar large velocity gradient solution showing 
the lowest kinetic temperature solution (100 K) for the cloud NGC 604-2 located in the giant 
HII region NGC 604. The values plotted for the ^^CO/^^CO J=l-0 and J=2-l ratios are the 
same as in (a); for the ^^CO J=3-2/J=2-l ratio, the plotted values are 1.07 ± 0.16, where 
the uncertainty includes the absolute calibration uncertainty. 

Fig. 3.— Large velocity gradient solutions with a kinetic temperature of 15 K for three 
different assumed values of the [^^CO]/[^^CO] abundance ratio. This hgure demonstrates 
the existence of solutions with different values of the abundance ratio but the same density, 
kinetic temperature, and ^^CO column density. The lines indicate the observed ila values 
for the line ratios of: ^^CO/^^CO J=l-0 (9.1,10.3, dashed line); ^^CO/^^CO J=2-l (6.8,7.8, 
dotted line); ^^CO J=3-2/J=2-l (0.63,0.75, solid line). The vertical solid line indicates the 
common density solution of 8300 cm“^. (a) p^CO]/p^CO]=30. (b) p^CO]/p^CO]=50. (c) 

p2cO]/[i3co]=70. 



Table 1. Log of JCMT Observations 





’^co and ’^CO J=2-l data 

’^CO J=3-2 data 

Cloud 

a(1950) 

(5(1950) 

Date'’ Time Time 

Date’’ Time 


m 

(O , 

(’2CO, m) (’’CO, m) 

(to) 


MC 1 

01:31:03.0 

30:23:56 

94/7/17'= 

18 

72 

95/8/1 

60 

MC 13 

01:31:10.3 

30:20:40 

94/7/16-17'= 

24 

62 

95/8/1 

60 

MC 19 

01:31:13.4 

30:23:52 

93/8/1 

24 

60 

95/8/1 

60 

MC 20 

01:31:10.9 

30:25:24 

93/8/1 

24 

72 

95/7/31'’ 

48 

MC 32 

01:30:51.8 

30:23:52 

93/7/30 

30 

42® 

95/8/1 

60 

NGC 604-2 

01:31:44.2 

30:31:27 

95/7/31 

18 

48 

95/7/31 

30 

NGC 604-4 

01:31:45.1 

30:30:58 

93/7/31 

18 

42 

95/8/1 

60 


^'Observations 1993 July 30 - August 1 obtained with the Canadian Acousto-Optical Spectrometer 
with 0.248262 MHz resolution; other observations obtained with the Dutch Autocorrelation 
Spectrometer (DAS) with 0.3125 MHz resolution. 

’’Observations obtained with DAS with 0.078125 MHz resolution. 

'’Telescope efficiency was 70% of normal. 

‘’Observations obtained with 0.3125 MHz resolution. 

^Observations obtained on 93/8/26 with 0.15625 MHz resolution. 







Table 2. CO Line Ratios in Molecular Clouds in M33 


Cloud 

I(i2C0 2-1) 

(K km s“^) 

AVfwhm 

(km s“^) 

i^CO J=2-l 

i^CO J=2-l 

12CO J=3-2 

i^CO J=2-l 

i^CO 4=1-0'^ 

i^CO J=l-0 

12±log(0/H)'^ 

HII region: 







MC I 

6.7 ±0.3 

8.0 

6.9 ± 1.0 

0.87 ±0.07 

9.2 ±2.0 

9.02 

MC 13 

6.8 ±0.4 

7.5 

6.2 ±0.7 

0.82 ±0.08 

10.1 ±0.6 

8.83 

MC 20 

13.6 ±0.2 

10.5 

7.4 ±0.6 

0.73 ±0.03 

9.0 ±0.4 

8.91 

NCC 604-2'= 

14.6 ±0.1 

11.0 

9.9 ±0.6 

1.07 ±0.03 

11.9±2.0 

8.51 

No HII region: 







MC 19 

8.2 ±0.4 

10.0 

6.3 ±0.4 

0.54 ±0.02 

11.3±0.8 

8.98 

MC 32 

12.2 ±0.3 

11.0 

8.2 ± 1.2 

0.70 ±0.02 

8.9 ±0.9 

8.83 

NCC 604-4 

11.8 ±0.4 

9.0 

6.4 ± 1.0 

0.50 ±0.04 

7.4 ±0.6 

8.51 


Uncertainties are uncertainties in the mean 

®'From Wilson & Walker 1994 

'^Metallicities interpolated from Vilchez et al. 1988 


This cloud is located in the giant HII region NGC 604 
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Table 3. Large Velocity Gradient Solutions for M33 Clouds 


Cloud 

pCO]/pCO] 

Tk 

(K) 

nH2 

(cm“^) 

iV(i2C0) 

(cm“^) 

NGC 604-2 

30 

200-300 

2 - 3 X 10^ 

20 X 10^^ 


50 

100-300 

1 - 3 X 10^ 

30- 70 X IQi^ 


70 

100-300 

1 - 3 X 10^ 

50 - 100 X 10^^ 

with HII regions 

30 

30-100 

2 - 6 X 10^ 

5-9 X 10^^ 

(MC 1, MC 13, 

50 

20-50 

2 - 7 X 10^ 

7- 20 X 10^^ 

MC 20) 

70 

15-30 

3 - 10 X 10^ 

7- 20 X 10^^ 

without HII regions 

30 

10-20 

5 - 30 X 10^ 

2-4 X 10^^ 

(MC 19, MC 32, 

50 

10 

10- 20 X 10^ 

3 X 10^^ 

NGC 604-4) 

70 

10 

10 X 10^ 

4 X 10^^ 

average of six clouds 

30 

15-30 

3 - 10 X 10^ 

2 - 6 X 10^^ 

(without NGC 604-2) 

50 

10-20 

4 - 30 X 10^ 

3 - 6 X 10^^ 


70 

10-20 

6 - 20 X 10^ 

4-6 X 10^^ 















































